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Abstract 

We propose that macroscopic objects built from negative-permeabihty metamaterials may ex- 
perience resonantly enhanced magnetic force in low- frequency magnetic fields. Resonant enhance- 
ment of the time-averaged force originates from magnetostatic surface resonances (MSR) which are 
analogous to the electrostatic resonances of negative-permittivity particles, well known as surface 
plasmon resonances in optics. We generalize the classical problem of MSR of a homogeneous ob- 
ject to include anisotropic metamaterials, and consider the most extreme case of anisotropy where 
the permeability is negative in one direction but positive in the others. It is shown that deeply 
subwavelength objects made of such indefinite (hyperbolic) media exhibit a pronounced magnetic 
dipole resonance that couples strongly to uniform or weakly inhomogeneous magnetic field and pro- 
vides strong enhancement of the magnetic force, enabling applications such as enhanced magnetic 
levitation. 



I. INTRODUCTION AND MOTIVATION 



Electromagnetic (EM) forces are used in a variety of technologies, ranging from motors, 
EM gun^ and magnetic levitation trains^ to micro- and nano-actuatord^ and optical tweez- 
xhe spatial scales where EM forces are important to mankind range from atomic to 
cosmic®, and the relevant EM frequencies cover the entire spectrum from statics to gamma- 
ray frequencies. 

For applications dealing with human scales (~ 1 m), high-intensity EM fields necessary 
for the generation of strong EM forces interfere with other devices and may be harmful 
to biological tissues. The severity of this problem is substantially reduced if the fields 
are predominantly magnetic, since virtually all biological substances and the majority of 
conventional materials are almost purely dielectric and thus transparent to magnetic fields. 
While at a finite frequency it is impossible to completely eliminate the electric fields from the 
picture, the ratio E/H can be strongly suppressed in the near field zone of a magnetic dipole 
oscillator, such as a conductor loop carrying alternating current (AC). For that reason, a 
large variety of EM force applications, including EM motors, maglev trains and magnetic 
brakes, have a strong preference for magnetic rather than electric forces. 

Natural media with strong magnetic response, although readily available, provide only 
a limited range of magnetic permeability (/i) values; negative permeability is particularly 
hard to find. Negative can be found in magnetized ferrites in narrow frequency bands 
associated with ferromagnetic!^ and gyromagneticP resonances. However, magnetic perme- 
ability typically has a large loss tangent in those bands, leading to high specific absorption 
rates. As we show in this paper, negative fi with low enough loss tangent enables resonant 
enhancement of the magnetic polarizability; therefore, it is important to find or construct 
low-density media with this property. 

Artificial media such as metamaterials have been demonstrated to provide negative fi with 
relatively low loss tangentd^^. Magnetic metamaterials consisting of metallic lines on thin 
dielectric substrates can be readily and inexpensively fabricated using the existing printed 
circuit board technologies (PCB). PCB metamaterials consisting of thin metal lines on low- 
density dielectric substrates have a substantial volume fraction of air and thus typically 
have much smaller density (of order 0.1 g/cm^) than ferrites and other naturally magnetic 
media. Their loss tangents can be reduced by using low-loss dielectric substrates and high- 
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conductivity metals. Additionally, the highly porous structure of PCB metamaterials allows 
for convection cooling, enabling higher peak intensities that do not cause material damage, 
which is impossible with homogeneous natural magnetics. The operational band of negative- 
fi PCB metamaterials can be easily adjusted to any frequency between roughly 1 MHz 
to about 100 GHz, even without conceptual changes in the PCB fabrication process, by 
adjustments to the unit cell size, as well as inductance and capacitance of the resonating 
metamaterial elements. 

Resonant enhancement of electric dipole polarizability by virtue of the so-called Surface 
Plasmon Resonance (SPR) is well-known in optic^^^^. Its closely related analog. Surface 
Phonon Polariton Resonance (SPPR) can be found in mid-infrarecP^'^. Surface plasmon 
resonances are electrostatic in nature and can be understood as eigenmodes, that is, source- 
free solutions, of the electrostatic Laplace equatiorpHSHSl 

VeV0 = 0. (1) 

When all materials involved have positive permittivity e > 0, a uniqueness theorem can be 
proven that precludes the possibility of non-trivial solutions to this equatioEpSIZI]_ However, 
when two or more homogeneous media with both negative and positive e are included, 
source-free solutions exist at certain, negative values of the permittivitj^^^**^. In binary 
(two-medium) composites, these values depend exclusively on the shape of the positive- 
negative e boundarj^^SHUl^ ]^-^J^^ l^j^g absolute size of the structure, since electrostatic 
problems do not have a spatial scale associated with the wavelength. 

For simple geometric shapes, resonant values of — e are not very large, of order unity: for 
a sphere, the dipole resonance occurs at Cres = —2, for a cylinder, it occurs at eres = — 1- 
The value of —eres can be increased by one or two orders of magnitude by introducing large 
aspect ratios and/or small gaps between the negative-e particles. Still, relative permittivity 
in the range —100 to —1 is readily found in natural media only at optical (UV through 
mid-IR) frequencies, which explains why surface plasmon resonances have so far been a 
subject of mostly optical studies, and attracted little attention from the point of view of 
RF and microwave applications. On the other hand, metamaterials can easily deliver those 
moderately negative values of e at any frequency of interest, including the microwave and 
lower-frequency regions. 

Magnetostatic fields in the absence of electric currents can be described by the magneto- 
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static potential which satisfies the equation V/iV</'*^™'' = 0, mathematically identical to 
the Laplace equation of electrostatics ([T]). This analogy is a particular case of EM duality, 
which can be stated as the exact equivalence of electric and magnetic field equations in 
the absence of free electric (and magnetic) charges. We can thus duplicate the electrostatic 
resonance phenomena using negative permeability media. This conclusion would also apply 
in the quasistatic, finite-frequency long as the retardation parameter D/X, where D 

is the characteristic spatial scale of the problem, is negligibly small. 

Existence of electric charges and their currents however creates a distinction between elec- 
trostatic and magnetostatic resonances. While the negative-e bands found in natural and 
artificial media can extend to arbitrarily low frequencies (recall the Drude dispersion of elec- 
tron fiuid), one cannot achieve negative /i(w) at zero frequency^, as the latter would imply 
negative energy density, which violates the conditions of thermodynamic stability. Nega- 
tive permeability always occurs in finite-bandwidth portions of the spectrum surrounded by 
positive-yU bands; in that regard, magnetostatic surface resonances (MSR) studied in this 
paper resemble quasi-electrostatic surface resonances in plasmon-polaritonic structure^^^^. 
Within the already narrow negative-/i band, we predict a series of even narrower resonances 
associated with magnetic surface eigenmodes of the macroscopic magnetizable body. 

In this paper, we consider metamaterial bodies that are much smaller than the free-space 
wavelength A of EM radiation at the operational frequency: D X. At the same time, 
we require that the unit cells of the metamaterial are much smaller than the whole body, 
a <^ D. This allows us to treat it as a homogeneous permeability object, and use metama- 
terial properties calculated from standard frequency-domain homogenization theories, both 
quasistatic!^ and full-wav^211_ xhe only currently known method for achieving negative-/! in 
metamaterials is the embedding of self-resonant metamaterial resonators, which thus need 
to be much smaller than D. The hierarchy a ^ D ^ A suggests that resonators must be at 
least two orders of magnitude smaller than A, which itself poses a challenging design issue. 
Fortunately, recent progress in microwave and RF metamaterials provides evidence that 
extremely subwavelength (A/a > 100) negative-/^ metamaterials can be manufacturecP^^. 

In studying magnetic surface resonances of metamaterial bodies, we cannot avoid dealing 
with the anisotropy of magnetic permeability. Magnetic response in metamaterials is created 
by resonant elements, which are more likely than not to respond only to one field polarization. 
Thus, from the design perspective anisotropic metamaterial^ are much easier to achieve 
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than their isotropic counterpart j^^ * ^^ l Additionally, magnetic resonators can be very thin, 
planar elements providing strong magnetic response only in the direction normal to their 
plane. A well-known example of such metamaterial elements are Split Ring Resonator^. 

Thin resonators can be stacked with high number density, providing stronger magnetic 
response for the select orientation of the magnetic field. The figure of merit for many applica- 
tions including superlens-based imagin^^^'^^'^ and wireless powei'^^'^is the inverse magnetic 
loss tangent at the frequency where Re/x(co') = — 1. Assuming the magnetic dispersion modeP 



/^(c^) = 1 - . ^ . X ^, (2) 



the aforementioned figure of merit is approximately 

Re/i 



Im/i 



IfQ, (3) 



where Q = Ures/i'^l) ^ 1 and F ^ 1. The magnetic oscillator strength F < 1 is propor- 
tional to the oscillator number density; thus, densely stacked thin resonators can reduce the 
loss tangent in metamaterials. 

By considering metamaterials with single-polarization response, we are naturally led to 
the consideration of magnetostatic (and by duality, electrostatic) resonances in indefinitd^ 
media. Electromagnetic waves in indefinite media have a hyperbolic dispersion relation!^, 
leading to a variety of interesting phenomena, including photonic van Hove singularities 
and extremely large photonic density of stateJ^. Since indefinite media have an unusually 
high density of states, applying the known eigenmode methods to them may be difficult, 
as discussed in the subsequent section. Already in driven EM problems, indefinite media 
can present numerical challenges, such as ill-conditioned stiffness matrices arising in the 
Finite Element Method. It is therefore both surprising and practically useful that indefinite- 
metamaterial bodies support isolated dipole resonances that couple strongly to uniform 
magnetic fields. In the final section of this paper, we show how these resonances can be 
utilized in magnetic levitation applications. 




FIG. 1. (color online). Schematic of the cylindrically-symmetric structure and the simulation 
domain used for calculations magnetostatic surface resonances. The PEC wire domain is included 



only in the magnetic force calculations described in Sec. Ill 



II. MAGNETOSTATIC SURFACE RESONANCES IN BODIES WITH NEGA- 
TIVE OR INDEFINITE MAGNETIC PERMEABILITY 



In exact analogy with electrostatic resonances, magnetostatic surface resonances can be 
obtained as source-free solutions of the equation 



V/xVc/.^'") = 0, 



(4) 



where c/)*-™-' (f) is the magnetostatic potential. In what follows, the superscript (m) is omitted. 
For the particular case of a homogeneous, possibly anisotropic, magnetic object with per- 
meability fj.^"^\co) occupying a finite domain ^l^n embedded in a homogeneous host medium 
with a permeability tensor (x^'^\uj), permeability can be expressed through the ^-function of 
the domain Q^: 

f,ix,oj) = + - e{x), (5) 

where by definition 6{x) equals 1 if x G Qm, and zero otherwise. 

For simplicity, we assume here that jj^'^^ = Hh is isotropic, and /i*^*"^ = diag(/ia;, /iy, /i^). In 
metamaterials, one can control all three principal values of individually; thus, the 

following Sturm-Liouville (eigenvalue) problem can be postulated: 



(6) 



where A = /i^ is the unknown eigenvalue and ^x,y are constants. From now on, we assume 
that the magnetic object is uniaxially anisotropic {fj.x = f'y), in which case Eq. ^ simphfies 
to 

/i;,V(l - e)V(l) + /i..V,,^V,.,0 = -A9,^a,0, (7) 
where Vxy = xdx + ydy. Alternatively, one may postulate an eigenvalue problem, 

/ift,V(l - e)V(t) + fizd,edz(p = -XVxyOVxy(j), (8) 

in which //j, = /i^ = A is treated as the eigenvalue. Finally, if the medium is isotropic, 
postulating fi^ = fiy = fiz = X results in the following eigenvalue problem: 

fihV{e - = AWV0. (9) 



Equations ( 6|[7p|[9 ) belong to the class of generalized eigenvalue problems (GEP), that 



is, they can be written as A(j) = XB(j), where A and B are linear operators; the particular 
case with 5 = 1 is known as the standard eigenvalue problem. Analytical properties of 
GEP can be improved by a rational transformation of the eigenvalue, A = {aX' + b)/ (A' + c), 
which leads to a new GEP A'(f) = X'B'(j), with the new linear operators A' = Ac — Bb and 
B' = aB — A. With a proper choice of a, b, c one can sometimes make the new operators A' 
and B' positive definite, and also confine all eigenvalues to a finite interval such as (0, 1). 
For instance, transformation A' = 1/(1 — fiz/f^h) converts Eq. Q to the canonical form 
postulated by Bergman et al.'^ 

VeV0 = A'V^. (10) 



The eigenvalues of the problem (10) are known^° to be real and confined on the interval 
(0, 1), that is, in the region(s) with Re[/i2(w)//U/i(w)] < 0. In what follows, we assume = 1- 
The eigenvalues can be found analytically only for several geometries that possess a high 
degree of symmetry, such as circular and elliptic cylinderd^H^ spj^ered^, hemisphere^ and so 



on. For most other shapes, the problem (10) must be solved numerically. However, once the 
eigenmodes of Eq. (10)) are found, analytically or numerically, the response of the structure 
with a complex permeability fi^ = l^'z + ^A^" to a uniform magnetic field Hq can always be 
written in the following form, no matter how complicated the shape of the structure: 

m = Mr) + E Ti^l^Mr). (11) 

„ X' -X'^ (0n, 0n) 



where the complex permeabihty enters as A' = 1/(1 — fiz/fJ'h), the functions 4>nif^ are the 
eigenmodes corresponding to the eigenvalues (each labeled with an integer subscript n), 
and 00 = ~HqZ is the magnetostatic potential of the uniform magnetic field polarized in the 
2;-direction. The scalar product in (11) is definecP^as {(f), ip) = J^6V(f)* ■ Vip- Additionally, 
as shown in Ref.'^, the total (magnetic) dipole moment of the inclusion Qm is always given 
by the expression 



m. = H,vY,j^, (12) 



no matter how complicated the shape of Qm is; general expressions for /„ in terms of the 
eigenmodes of each resonance can be found in Ref."^. Thus, calculation of a magnetostatic 
response of a given body can be reduced to the calculation (analytical or numerical) of the set 
of numbers and /„ that depend only upon the shape of Qm and characterize its magnetic 
polarizability completely. The numbers A'„, different for each unique shape, determine the 
values of the magnetic permeability fiz at which magnetostatic resonances occur, and /„ 
determines the strength of that resonance. 

A general observation regarding the shape of the magnetic response as a function of 



permeability of the body can be made from the equation (12). As a function of the complex 
parameter A' = 1/(1 — fiz) = A^ + iX'^, magnetic moment is a sum of simple-pole terms. 
For any passive medium, Im/i^ and consequently ImA' = A^ is positive; as we pointed out, 
the eigenvalues A'„ are always reaP^. If we assume a constant value of A'j > and plot the 



real and imaginary parts of one resonant term from Eq. 12, as a function of the real 

parameter A^, we trivially obtain 

1 A' — A' 

^'a^ = {K-h' + iKr ^^^^ 

Comparing these expressions with the well-known shape of resonant response functions in 
the frequency domain, i.e., 

u"^ — u"^ — i'ju (w^ — tuy + (7cj)2 ' 

— — i'jco (w^ — u;)^ + (70;)^ ' ^ ^ 

we must conclude that 77Z^ clS db function of A', looks like the familiar Lorentz-shaped res- 
onance. As a function of /i'^, however, the magnetic polarizability deviates from a perfect 
Lorentzian shape, because A'^ = (1 — /i'^)/[(l — /x'^)^ + (/x")^] is not a linear function of fi'^. 
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For an arbitrary-shape particle Qm, the aforementioned eigenvalue problem must be 



solved numerically. The GEP (10) describing the isotropic medium case was solved in 
j^gfgEHIS] using the Finite Element Method. We have generalized the method of Ref.'^ 
to include the anisotropic medium equations ( 6||7||8 ). Numerically, the obtained eigenvalue 



problem produces reasonable results whenever the tensor fi^^\u) is negative-definite, i.e. 
when jjix = jJLy < Q and /i^ < 0; in that regime, both Eqns. (7||8) returned a sequence of 



well-separated eigenmodes. However, in the most interesting, indefinite medium regimes, 
including fi^ = jjy < 0, jj^ > and fi^ = fiy > 0,Hz < 0, both Eqns. (7]|8) produce an 



extremely dense set of eigenmodes, all of which have a very small dipole moment, and cor- 
respondingly negligible coupling to uniform magnetic field. Decreasing the mesh size in the 
FEM discretization results in even denser sets of eigenmodes; thus, no /^-convergence can be 
observed. 

This numerical problem stems from the well-known fact that indefinite media support 
waves with a hyperbolic dispersion relatiorP^. A lossless indefinite medium thus allows wave 
propagation with arbitrarily large wavenumberd^ in practice, there is still a maximum 
(cut off) wavenumber kmax, which is related to metamaterial granularity and loss, but is 
typically much larger than the free-space wavenumber k^. In numerical FEM solutions, 
kmax is set by the mesh size h; shrinking the latter leads to an even larger number of 
eigenmodes corresponding to high /c-waves propagating in the indefinite medium. Since no 
/i-convergence can be obtained, there is no straightforward generalization of the method 
from Refs.'I^^ to the indefinite medium case. 

On the other hand, it is known that at least simple shape objects filled with indefinite 
medium support well-defined dipole resonances, which couple strongly to the uniform field. 
For example, the magnetic polarizability of a homogeneous sphere with relative permeability 
fi = diag(/ia;, /iy/iz) can be found by applying EM duality to the electric polarizability derived 
by SihvolaP: 

where V is the volume of the sphere. Therefore, even with fi^ = fiy > 0, there exists a 
strong, distinct and isolated dipole resonance at /i^ = —2, situated deeply in the indefinite 
medium regime. Analytical solutions of the surface resonance problem are known only for a 
handful of simple shapes, including circular and elliptic cylinders and spheres. To the best of 
our knowledge, the sphere is the only essentially three-dimensional shape of an anisotropic 
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FIG. 2. (color online). Real (a) and imaginary (b) part of the volume- normalized magnetic 
polarizability ai™^ of a metamaterial disc, versus the value of Re/x^ in the disc, for various shape 
aspect ratios {H/R, indicated by numbers in the legend). Peaks of Imai'J*^ correspond to MSRs 
of the disc with non-vanishing dipolar strength. The disc is homogeneous with permeability fi = 
diag(l + O.li, 1 + O.li, Re/i^ + O.li). 

medium object for which analytical solutions were obtained. 

In order to find the eigenmodes of equations 
metamaterial body, we resort to solving a driven magnetostatic problem. The numerical 
model includes a homogeneous (effective) medium model of a metamaterial body and a 
finite-size simulation domain of diameter substantially larger than the metamaterial object, 
but still much smaller than the free-space wavelength. A constant, uniform magnetic field 
Hq = HqCz is applied, which corresponds to magnetostatic potential 0o = ~Hoz; here, 
-f^o = 1 A/m is an arbitrary intensity of the applied magnetic field. The metamaterial 
permeability is then scanned through the range of values where magnetostatic resonances 
are expected. For brevity, we describe the solutions only to the problem ([T]), in which /i^ 
is treated as the unknown (eigenvalue) and /i^. = /i^ = 1; the other eigenvalue problems 
can be treated similarly. To simulate the response in the z-direction, we choose the exterior 
boundary of the simulation domain to be a circular cylinder with its axis of revolution 
pointing in the z-direction. Dirichlet boundary conditions = and = —HoHdom are 
applied on the bottom {z = 0) and top {z = Hdom) plates; the side walls are modeled as 
Neumann boundaries (n ■ V0 = 0), which implies that the magnetic field is parallel to those 
boundaries. 



6|[7p) for an arbitrary-shape indefinite 
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For the composite structure consisting of the metamaterial body and an homogeneous 
isotropic medium surrounding the metamaterial, one can introduce the relative effective 
permeability in the direction of the applied magnetic field: 

/xf = , ^ ^ / n . BdA, (18) 

where Adom = T^R^dom the area of the bottom boundary of the cylindrical domain, and 
integration is carried over the latter. The simulation domain for this type of calculation 
is depicted in Fig. [Tj This is the dual, magnetic version of the "capacitor" definition of 
effective permittivitj^^^. In what follows, the metamaterial is assumed to be situated in free 
space = !)• A finite loss characterized by Im/x'^™'^ > is assumed in the simulations; 
in the indefinite medium case, this turns out to be necessary not only for physical realism, 
but also to provide numerical stabilization of the resulting linear algebra problem. Without 
the loss, the stiffness matrix of the FEM discretization can be ill-conditioned, i.e. close to a 
singular matrix. 

In the magnetostatic limit, one can readily relate the effective magnetic permeability 



as defined by Eq. 18, to the magnetic polarizability of the magnetizable body inside the 
simulation domain. Integration of the local magnetization density, M = B/fiQ — H, over the 
simulation domain Q, yields the total magnetic moment of the body, since it is situated in 
free space where M = 0: 

[ MdV = I MdV = m. (19) 

One can then introduce effective susceptibility averaged over the domain Q in accordance 
with the following definition: 

~ HoV ~ HoV ~ V ' ^ ^ 

where we have also introduced the magnetic polarizability, aiT"* = mz/Ho, of the magneti- 
zable body represented by the domain Q^n ^ ^- In terms of dimensionless quantities - the 
filling fraction = Vm/V and volume-normalized polarizability ai?-* = ai^V^' where Vm 



is the volume of domain Qm - Eq- 20 can also be stated as xl^ = fmC^zf' ■ 



It was shown by Urzhumov et al. using the terminology of electrostatic^ that, in the 
quasistatic limit, the quantities defined by Eq. ( 18p0 ) are related by the familiar equation 



/^f - 1 = Xf , (21) 
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FIG. 3. (color online). Magnetostatic resonances of a homogeneous disc from Fig. [2] and aspect 
ratio Hciisc/ Rdisc = 0.5. Cross-sections show magnetostatic potential cj)^'^ = cp — (pQ, calculated 
from Laplace equation with excitation potential 0o = —z, corresponding to uniform magnetic field 
of unit magnitude in the z-direction. (a) Potential distribution at Re/x^ = —0.49. (b) Potential 
distribution at Re/x^ = —1.68. 



and consequently, ^'f — 1 = fmC^zz'- Thus, the calculation of fi'f — 1 by virtue of Eq. 



cff 



also gives the magnetic polarizability of the body, which are directly proportional to each 
other. 

The resonances are identified as Lorentz resonances in Re/Xeff and IniyUefr plotted as a 
function of Re/i^ (see Fig. |2]). The metamaterial permeability is assumed to be /i = diag(l + 
O.li, 1 + O.li, Re/ix + O.li), with variable Re/x^. The finite metamaterial loss lets us locate and 
characterize the resonances that couple to uniform magnetic field by plotting Im^^s versus 
Re/i^; such a curve features a series of well-defined absorption bands. The abscissa of each 
absorption peak defines the value of Re/i2(a;) (and, given the metamaterial dispersion curve 
fiz{^), also the frequency) at which the resonance occurs; the ordinate measures the strength 
of coupling of the resonance to the uniform magnetic field. Fig. |2] shows the results for a disc- 
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shaped metamaterial body shown schematically in Fig. [TJ its cylindrical symmetry allows us 
to speed up the calculations by using the axisymmetric formulation of the magnetostatics 
problem in COMSOL Multiphysic^. 

In the range —2 < Kefiz < 0, we observe one or two magnetic-dipole resonances of the 
disc, which occur at different values of Re/Xz depending on the aspect ratio Hdisc/ Rdisc, 
where Hdisc and Rdisc are the height and radius of the metamaterial disc, respectively. In 
the entire range — oo < Re/i^ < 0, there is an infinite series of such resonances occurring 
at progressively higher values of — Re/i^; however only the lowest two resonances are strong 
enough to have practical applications. In particular, we do not find strong resonances 
associated with the edge of the disc; in general, sharp corners tend to produce resonanceP^ 
with fields confined around the tip of the corner. In the case of a disc-shaped object, which 
has only edges and no sharp corners, such eigenmodes do not contribute significantly to the 
quasistatic response. 

For the aspect ratio Hdisd Rdisc = 0.5, the first and the strongest resonance is obtained 
at Re/i^ ~ —0.5, and the second resonance — at Re/i^ ~ —1.65. The magnetic potential 
distributions corresponding to these resonances are shown in Fig. [3j The implications of 
these resonances to magnetic force enhancements are discussed in the next section. 



III. MAGNETIC LEVITATION OF DEEPLY SUBWAVELENGTH METAMATE- 
RIAL OBJECTS 



The potential energy of a magnetic body possessing a permanent magnetic dipole moment 
m, magnetic quadrupole moment Q^"^^ and other multipole moments can be written aP^l 



U{x) 



(22) 



the magnetic field is evaluated at the geometric center of the body. This multipolar expansion 
is valid whenever the body is substantially smaller than the spatial scale of magnetic field 
inhomogeneity link = |jv^||q|- Using the principle of virtual displacement, one derives the 
force 



dU 

dxi 



OH, 



(m) 



OH, 



dxj 6 ■^'^ dxjdx, 



+ 



(23) 



For bodies much smaller than Imh-, the first, dipolar term is typically the dominant con- 



tribution both in the potential energy (22) and in the magnetostatic force (23), and higher 
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FIG. 4. (color online). Magnetic field {\H\) distribution on the surface of a metamaterial disc with 
permeability n = diag(l + O.li, 1 + 0.1i,Kefiz + O.li). Disc aspect ratio Hdisc/ Rdisc = 0.5. (a) 
Re/x^ = -0.5, (b) Re/i^ = -1.65. 

order multipoles can be neglected. Assuming the geometry shown in Fig. [4], that is, a disc 
of radius Rdisc over a circular coil or radius Rcoih the scale Imh is on the order of Rcoii- 
For Rdisc ~ Rcoih quadrupolar contributions can be noticeable; to simplify our analysis, we 
assume Hdisc < Rdisc ^ Rcoih which allows us to focus on the lowest-order term in the 



expansion (23), which is proportional to the magnetic dipole moment. 



When the magnetic body has no permanent and only induced magnetic moment, one 



introduces the magnetic dipole polarizability tensor aj™'' such that m 



For 



induced moments, the potential energy expression (22 ) needs to be corrected by an additional 



factor of 1/2, however the force expression (23) is still correct, and it reads, in terms of the 
polarizability tensor, 

where the dots represent higher-order terms that are small relative to the dominant, dipole 
term for bodies of diameter <^ Unh- Although our main interest is to anisotropic bodies, we 
should remind the reader that for isotropic bodies (i.e., belonging to an isometric symmetry 



class and filled with isotropic permeability), magnetic polarizability a- 



(m) 
jk 



a 



Ma. 



fc IS a 



spherical tensor; in that case, Eq. (24) simplifies to the familiar expression F = ^a^"^^'V{H'^). 
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Note that this force vanishes in a uniform magnetic field, regardless of the properties of the 
magnetizable object. 

Before we consider magnetic levitation of metamaterial bodies, let us make a crude es- 
timate of the current in the coil required to levitate an object of an arbitrary, unoptimized 
shape, filled with isotropic, positive permeability. Magnetic polarizability of a sphere of 



volume V and isotropic relative permeability /i^ is, according to Eq. (17) 



= /ioV^^^^\^ ~ fioV, (25) 

where in the latter estimate we assume fir > 1- Magnetic induction on the axis of the 
coil carrying current I is = fioR'^^^iI / {2d^) , where d is the distance from the coil center. 
Assuming d ~ Rcou, we estimate the B-field gradient to be ~ B/Rcou, which gives the 
magnetic force density 

in SI units. Dividing this by the gravitational force density Fg/V = pg, where 10 m/s^, 
we obtain the figure of merit for magnetic levitation: 

Thus, levitation of an object of density 1 g/cm^ and size 1 cm or smaller using a single-turn 
coil of radius 1 cm requires a current of order 100 A, well beyond the practically achievable 
limits. It is therefore important to find a way to enhance the magnetic force by optimizing 
the properties of the levitating object. 

From the previous section we already know that tuning the permeability of an object 
to one of its magnetic surface resonances results in resonant enhancement of the magnetic 
polarizability, and consequently, one expects a similar increase in the magnetic force. In 
a finite-loss medium, the enhancement is proportional to the factor 1/lmfi'^^, which arises 



from the denominator of Eq. (25) in which Kefi^^ is matched to the resonance condition. 
Practically, this means that a resonantly-enhanced levitation system needs a smaller current 
magnitude — smaller by a factor of |Im/i'^'^|^^/^ — to achieve the same levitating force that 
would exist in off resonance conditions. Alternatively, with a given, attainable current mag- 
nitude, one could levitate an object heavier by the factor of 1/lmfi'^^. Loss tangents of order 
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0.1 are attainabl^^^'^in negative-permeability metamaterials at 10 MHz; thus metamaterials 
could potentially increase magnetic levitation strength by an order of magnitude. 

To quantify this phenomenon more robustly, we carry out full-wave EM simulations based 
on the Finite Element Method. We solve the Helmholtz equation for the vector potential 
A, which is linked to the magnetic induction in the usual fashion, B = V x A. Again, we 
assume a disc-shaped metamaterial object of radius Rdisc and height Hdisc levitating above 
a circular coil of radius Rcou and wire radius r^ire = 0. 01 Rcou- The simulation domain is 
shown schematically in Fig. [Tj When the disc is situated coaxially with the coil, we can take 
advantage of the cylindrical symmetry of the problem and use the axisymmetric version of 
the EM solver as provided in COMSOL Multiphysics^^. The vertical distance between the 
geometric center of the coil and the center of the disc is denoted d. 

The coil is driven by a fixed AC current source, which delivers the current /q = 250 A at 
frequency 10 MHz, corresponding to the free space wavelength Aq ~ 30 m. The metamaterial 
disc is modeled as a homogeneous body with permeability fi = diag(l + O.li, 1 + O.li, Re/iz + 
O.li), where Re/i^ can vary. The total force acting in the vertical direction (along the axis 
of revolution) is calculated by the software using surface integration of the Maxwell stress 
tensor; more details on this formalism can be found e.g. in Refs.'^^'^. Neglecting the already 
vanishingly small ratio Rdisc/ \i there are three dimensionless parameters describing this 
configuration, in addition to the disc permeability: the disc aspect ratio Hdisc/ Rdisc, disc 
size relative to the coil Rdisc/ RcoU, and finally the coil-disc distance relative to the coil size, 
d/ Rcoii- 

When the ratio Rdisc/ RcoU is small, the effect of higher-order multipole resonances can 
be neglected, which simplifies the picture and allows us to make a clear identification of 
the dipolar MSRs. For that reason, we assume Rdisc/ RcoU = 0.2 here. Note that both 
the magnetic force associated with dipole polarizability and the gravity force scale linearly 
with the volume of levitating body; therefore, choosing to work with smaller bodies neither 



improves nor decreases the levitation figure of merit (27). 

The dependence of the accurately computed cycle-averaged total EM force on the re- 
maining parameters, d/Rcou, Hdisc/ Rdisc-, and /i^, is reported in Fig. |5} First, we choose the 
optimum disc position relative to the coil by scanning d/Rcou- A well-defined optimum is 
situated aX d = OARcou for our choice of parameters, as seen from Fig. [sj^a). The existence 
of this optimum is as follows: at sufficiently large distances from the coil, the magnetic field 
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on the axis falls off as however, the force is also proportional to the field gradient. 

Due to symmetry, B-field gradient vanishes in the plane of the coil (at d = 0). Therefore, a 
location exists in the vicinity of the coil where \d{B'^)/dz\ has a maximum. 

From Fig. we observe that, as a function of Re/i^, the force resonance always occurs 
at a fixed value of Re/i^, regardless of the distance d. This proves that this resonance is 
entirely due to the properties of the metamaterial disc, and has nothing to do with the 
magnetic field configuration around it. In Fig. [5|b), we fix the coil-disc distance d at its 
optimum value, and study the position of the force resonance as a function of the disc shape 
aspect ratio, Hdisc/Rdisc- 

Comparing the force resonances seen in Fig. [5|b) to the resonances of magnetic polar- 
izability (i.e. — 1) of the metamaterial disc in Fig. [2| we observe a striking similar- 
ity. Both figures reveal a Lorentz-shape resonance, which is positioned at Re/i^ = —1 for 
Hdisc/Rdisc = 1, Re/i^ = -0.5 for Hdisc/Rdisc = 0.5, and Re/x^ ^ for Hdisc/Rdisc < 1- The 
resonance positions are in precise agreement between in Fig. ^h) and Fig. |2} This enables us 
to identify the strongest force resonance in Fig. ^h) with the strongest dipole MSR visible 
in Fig. [2] for the aspect ratio Hdisc/Rdisc = 0.5 it is also visualized in Fig. [3|a). 



IV. CONCLUSIONS 

In conclusion, we have theoretically demonstrated the possibility of resonant force en- 
hancement in metamaterial bodies with negative permeability. Force enhancement is linked 
to the magnetostatic surface resonances (MSR) which are analogous to localized surface 
plasmon resonances (SPR) in negative-permittivity particles. Resonances of the magnetic 
force are approximately Lorentz-shaped as a function of the magnetic permeability of the 
metamaterial object. Consequently, the magnetic force induced in negative-/i objects can be 
either repulsive or attractive, depending on the object shape and its magnetic permeability. 
The enhancement factor with respect to the positive-permeability case is given roughly by 
1/Im/i''^, where /i'^^ is evaluated at the resonance condition Kefi'^^ = fires- With realistic 
metamaterial loss tangents of 0.1, this phenomenon could allow magnetic levitation systems 
to increase the mass of levitating objects by one order of magnitude while using the same 
current magnitude. 
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Re Re 

(a) (b) 

FIG. 5. (color online). Cycle-averaged electromagnetic force in the z-direction acting on a meta- 
material disc with permeability ^ = diag(l + 0.1?, 1 + 0.1i,Re//2 + O.li), versus the value of 
Re//^ in the disc. Vertical force is normalized to the gravitational force Fg = VdiscPg, where 
y^disc = R%sc^ disc-, assumiug p = I g/cm^ and 5 = 10 m/s^. Positive force implies repulsion from 
the coil (lifting force). The current in the single-turn coil is 250 A; coil radius Rcoil = 2 cm and 
disc radius Rdisc = 0.2Rcoii = 4 mm. (a) Force variation as a function of the distance d (varied in 
the range 2-20 mm) between the geometric centers of the disc and the coil, for the disc aspect ratio 
Hdiscl Rdisc = 0.5. (b) Force variation as a function of the disc aspect ratio Hdisc/ Rdisc = 0.1 — 1, 
for d = OARcoii = 8 mm, which corresponds to the maximum force enhancement in (a). 
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